基于 虚拟 天 文 馆 对 全 天 相机 姿态 的 精确 测量 
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摘要 全 天 相机 已 广泛 地 使 用 在 气象 、 天 文 等 领域 . 在 监测 云 量 、 夜 天 光 、 流 星 等 应 用 中 , 均 需 要 准确 掌握 相 
机 的 姿态 参数 ， 特 别 是 在 安装 和 运 维 中 造成 的 偏差 . 本 研究 旨 在 开发 一 套 基 于 虚拟 天 文 馆 的 像 场 测 量 方法 , 通 
过 比较 参考 恒星 的 虚拟 坐标 及 其 像 点 的 实测 坐标 , 高 精度 测量 相机 的 姿态 参数 . 使 用 中 国 科 学 院 云 南天 文 合 安 
ha 置 于 四 川 无 名 山 观 测 站 的 全 天 相机 , 分 析 于 2016 至 2017 年 获得 的 4 组 全 天 图 像 ， 采集 了 容量 分 别 为 ~10、~50、 
pas ~200 的 恒星 样本 , 通过 Stellarium 虚拟 天 文 馆 获得 了 恒星 的 地 平 坐 标 , 对 全 天 相机 像 场 的 天 顶 位 置 、 测 者 子午 线 
N 方向 等 基本 参数 进行 了 精确 的 测量 . 主要 结果 为 : (1) 开 发 了 一 种 基于 虚拟 天 文 馆 计 算 恒 星 地 平 坐标 功能 的 测量 
O Fk, 使 姿态 测量 建立 在 容量 较 大 的 参考 恒星 上 . 该 方法 准确 度 高 ， 对 设备 配置 的 依赖 性 低 ， 具 有 较 强 的 可 移植 
O E (2) 采用 网 格 法 对 图 像 天 顶 的 定位 达到 了 亚 像素 的 精度 ， 有望 满足 高 定位 精度 监测 的 需要 . (3) 全 天 相机 姿态 精 
O ” 度 的 首要 指标 是 光 轴 的 坚 直 性 , 光 轴 偏离 天 顶 会 对 投影 的 轴 对 称 性 造成 不 可 忽略 的 影响 . (4) 维 护 作 业 对 圆 形 像 
场 几 何 参数 的 改变 甚 微 , 但 可 能 会 改变 指向 和 旋转 角 , 并 需要 重新 测量 . (5) 生 成 了 符合 制图 惯例 的 可 视 化 产品 ， 
为 监测 数据 的 深入 分 析 英 定 了 基础 . 
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Accurate measurement of all-sky camera's attitude 
based on virtual planetarium 


LI Xiaobo', LIU Yu', CUI Chenzhou’, SONG Tengfei', WANG Jingxing' , FAN Jiankang' 


1. Yunnan Observatories, Chinese Academy of Sciences, Kunming 650011, China 

2. National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China 

Abstract All-sky cameras have been widely used in meteorology, astronomy and other fields. In the monitoring of 
cloud cover, night skylight, meteors and other applications, it is necessary to accurately measure cameras’ attitude, 
especially the deviations caused during installation and maintenance. The purpose of this study is to develop a method 
to measure camera's attitude parameters with high precision, by comparing the polar coordinates of reference stars' 
image points and the stars' horizon coordinates computed by a virtual planetarium, and generate visualization products 
with errors corrected. Using the all-sky camera developed by the Yunnan Observatories of the Chinese Academy of 
Sciences and deployed in Sichuan's Wumingshan observation station, we analyze four sets of all-sky images obtained in 
October 2016 and June 2017, select star samples with capacities of ~10, ~50, and ~200, obtain their coordinates through 
the Stellarium planetarium, and accurately acquire the attitude parameters such as the zenith position and rotation angle 
of the image field. Main results are: (1) We developed a method that is based on virtual planetarium's function to 
compute stellar coordinates, and guaranteed the accuracy of the measurement with a large capacity of reference star 
samples. The method does not rely on the camera’s specs and therefore is transplantable. (2) According to the camera's 
attitude on 2016-10-20, its image field's radius and center coordinates are 0.4289 and (0.7590, 0.5048) respectively (in 
the unit length of the image's Y direction, the same below); the zenith's coordinates are: (0.75118+40.00020, 
0.49572+0.00022); the angle between the zenith and the optic axis is 2.36°; the rotation angle is 152.71+1.58. (3) 
Maintenance operations have little (<0.28%) influence upon the geometric parameters of the round-shaped image field, 
but might cause the optical axis to deviate from the zenith, which will have an non-negligible effect on the axial 
symmetry of the camera's projection. Each maintenance will also change the camera's rotation angle which needs to be 
re-determined afterwards. (4) Produced visualization products conforming to industry norms, thus cements data’s 
reliability for further investigations. 


Keywords all-sky camera; virtual planetarium; field of view; zenith; monitoring of cloud cover; 
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1 引言 


全 天 相机 是 利用 朝向 天 顶 的 鱼 眼 镜头 获取 全 天 空 图 像 的 照相 机 .全 天 相机 已 获得 广泛 的 应 用 ， 


领域 . 
蓝天 ， 然 后 进 


(1) 地 基 遥 感 全 天 空 云 量 自动 化 监测 023. 在 数字 图 像 处 理 软件 中 ， 利 
步 对 云 进行 分 类 和 计 
进行 判断 9 其 他 计量 云 量 的 方法 有 


虚拟 天 空 "”、 比 较 夜 间 实 测 和 虚拟 星 
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主要 在 以 下 
书 经 验 或 理论 确定 的 阔 值 判别 云 或 


量 ， 大 多 数 算法 中 都 是 基于 瑞 利 散 射 或 饱和 度 对 从 三 原色 计算 出 的 菜 指 标 
: 比较 天 空 辐射 亮度 、 比 较 两 个 或 多 个 窄 波 段 的 辐射 值 、 比 较 白 昼 实 测 和 
EA. HAAR!) | ERED) 等 . 


对 于 注 云 、 临 近 太阳 的 天 区 、 特 殊 天 


气 情况 、 夜 间 云 量 、 


临近 地 平 线 的 云 量 的 判断 仍 是 难点 . 对 云 高 和 云 状 的 


动 监测 仪器 仍 需 进一步 实验 研究 . (2) 


监测 天 文 址 点 的 夜 天 光 和 云 量 "'" 1 在 我 国 的 天 文 界 , 全 天 云 量 相机 已 应 用 
海 观 测 站 中 、 西 部 天 文选 址 等 . (3) 用 于 监测 流星 的 全 天 相机 网 络 "%. (4) 对 舰 船 进行 天 文 导 航 定位 


于 丽江 天 文 观 测 站 0 、SONG 项 目 青 


方向 . 


洁 等 作业 也 会 改变 其 姿态 . 


所 以 , 在 处 理 全 天 相机 图 像 的 工作 


在 以 上 各 方面 的 应 用 中 , 均 需 要 准 
图 像 中 的 实际 方向 和 天 顶 位 置 ， 以 确定 特 包 


理想 情况 下 , 全 天 相机 的 光 轴 指向 天 顶 , 视 场 中 心 与 天 顶 
的 全 天 相机 不 能 达到 这 个 标准 . 首先 ， 
像 安装 大 型 天 文 观测 设备 那样 对 方位 进行 准确 的 校准 , 难 
无 人 值守 的 , 恶劣 天 气 、 时 生动 物 等 外 力 都 有 可 能 使 相机 


确 掌 握 像 场 的 参数 . 较 高 标准 是 : 准确 获得 物 一 像 的 投影 关系 ,能够 将 图 
象 中 的 特征 点 还 原 为 天 球 上 的 坐标 , 这 对 于 监测 和 定位 航空 器 、 流 星 等 小 尺度 物体 至 关 重 要 . 较 低 标准 是 : 掌握 


相机 一 般 安 置 于 高 山 、 


F 物 的 方位 . 对 于 任何 用 途 的 全 天 相机 图 像 ， 起码 要 知道 其 东南 西北 


EA; 图 像 的 纵 轴 与 测 者 子午 线 重合 . 但 是 ,实际 
沙漠 、 建 筑 物 顶 端 等 ， 
免 产 生 不 能 忽略 的 安装 偏差 . 其 次 , 全 天 相机 大 多 是 
局 离 其 初始 安装 姿态 . 再 


由 于 条 件 所 限 ， 施 工时 无 法 


次 , 对 相机 的 维修 、 改 造 、 清 


中 ,就 需要 依据 参照 物 对 图 像 进行 矫正 ， 获 得 其 旋转 


和 天 顶 位 置 等 基本 


参数 . 不 仅 需 要 对 新 安装 的 相机 初 光 获 得 的 图 像 进 行 精细 测量 ,在 日 常 运行 中 也 要 不 定期 再 次 测量 并 更 新 参数 . 
如 此 方 能 生成 可 靠 的 可 视 化 产品 , 为 监测 工作 提供 有 效 的 服务 . 另外 ,准确 测量 像 场 也 是 矫正 镜头 畸变 、 进 行 坐 


标 变换 、 进 而 实现 上 述 “ 较 高 标准 ”的 起 点 . 


恒星 是 全 天 相机 视 场 中 最 可 靠 的 参照 物 ， 而 虚拟 天 文 馆 大 大 简化 了 对 星 表 的 检索 和 计算 工作 , 在 图 形 化 界 


直观 的 用 户 界面 ， 往往 可 以 离线 工作 
的 实际 条 件 . 


较 特 征 恒 星 的 地 平 坐标 及 其 像 点 在 像 引 


外 中 可 以 快速 完成 检索 恒星 和 测量 坐标 的 操作 ,正好 可 以 用 


有 深入 开发 该 方法 . (2) 基 于 一 台 


及 其 与 像 场 中 心 的 距离 ,达到 亚 像素 上 


NG; 


次 重大 维护 之 前 和 之 后 分 别 涡 


i 


量 上 述 参量 ， 了解 维 护 工 作对 相机 姿态 的 影 


的 可 视 化 产品 . 下 面 介绍 本 研究 的 研究 资料 、 数 据 压 缩 方 法 和 研究 结果 


于 对 全 天 相机 姿态 的 测量 . 同时 ， 
， 对 计算 机 配置 要 求 相对 较 低 ， 相 比 于 虚拟 天 文 台 20， 更 符合 野外 观测 站 


虚拟 天 文 馆 具 有 


基于 上 述 需 求 和 条 件 ， 本 研究 拟 达 到 四 方面 的 目标 : (1) 开 发 一 套 基于 虚拟 天 文 馆 的 像 场 测量 方法 , 通过 比 
PH ERIKA, 高 精度 测量 相机 的 姿态 ， 并 探索 在 上 自动化、 可 移植 化 等 方 
型 全 天 相机 的 实际 监测 数据 ,通过 像 场 测量 其 姿态 : 确定 图 像 中 天 项 的 位 置 ， 
测量 测 者 子午 线 的 方位 ; 测量 相机 像 场 的 基本 几何 参数 . (3) 在 某 


响 . (4) 矫 正 姿态 偏差 并 生成 符合 规范 
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图 1 无 名 山 全 天 相机 于 2017 年 6 月 5 日 日 初时 拍摄 的 一 帧 图 像 
Fig. 1 An image taken by theWumingshan all-sky camera at the dawn on June 5, 2017. 


2 研究 资料 


本 研究 使 用 的 全 天 相机 架设 于 四 川 稻城 县 的 无 名 山 二 号 观测 站 . 无 名 山 是 中 国 科学 院 云南 天 文 台 选 址 与 日 冕 观 
测 团 组 在 西部 天 文选 址 科 考 工作 中 发 抉 出 来 的 优秀 天 文 址 点 > 于. 该 相机 由 光学 、 采 集 、 控 制 、 和 保护 系统 
组 成 . 其 中 ,光学 系统 基于 一 枚 适马 (Sigma) 牌 焦距 4.5mm、 光 圈 2.8 的 鱼 眼镜 头 ， 视 野 略 大 于 180°; 采集 系统 基 
于 一 台 佳 能 (Canon) 牌 EOS-60D 单 镜头 反光 相机 ; 控制 软件 和 数据 库 安装 于 一 台 研 祥 牌 工业 控制 计算 机 ,根据 日 
月 的 高 度 、 月 相 、 天 光亮 度 调整 感光 度 和 曝光 时 间 ; 光学 和 采集 系统 被 封装 于 一 个 能 抵御 恶劣 天 气 的 机 箱 中 , 镜 
头 透 过 透明 半球 单 指 向 天 顶 ， 形 成 圆 形 的 像 场 (图 1). 图 2 展示 了 位 于 机 柜 顶 部 的 天 窗 的 基本 结构 ， 该 半球 尝 半 

径 为 10 cm， 镜 头 前 端 高 出 半球 球 心 5 cm, 视野 达到 地 平 线 以 下 . 在 半球 单 的 周围 安装 了 用 于 除 箱 露 的 电热 装置 . 
自 2015 年 11 月 初 光 以 来 , 该 设备 对 无 名 山 的 天 空 进行 了 连续 稳定 的 监测 . 观测 人 员 对 半球 罩 进 行 例 行 清洁 , WF 
发 人 员 不 定期 对 相机 进行 检修 、 更 换 零 件 、 升 级 部 件 等 操作 . 


T 
Coy 


ø 


图 2 无 名 山 二 号 点 全 天 相机 的 设计 图 ( 左 ) 和 实际 外 观 ( 右 ) 
Fig. 2 The design (left) and actual outlook (right) of the all-sky camera deployed in Mt. Wumingshan No.2 observation station. 


Stellarium 


Pole ph Fabien Cháreau 等 人 开发 的 一 种 开源 虚拟 天 象 模拟 软件 . 基本 星 表 包 含 60 万 余 颗 恒星 ， 附 
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加 星 表 包含 177 万 余 颗 恒星 ， 
基于 无 名 山 二 号 观测 


Z] 


获得 该 时 该 地 的 星 图 ， 


站 的 经 纬度 (纬度 :29°6'24.59"”、 


(HIP21502) 的 详细 参数 , 包括 在 指定 观测 条 


H, PE. 


拟 分 析 无 名 山 全 天 相机 于 2016 年 10 月 和 2017 年 06 月 获 
好 的 气象 条 件 下 拍摄 ,用 于 分 析 星 场 ; cee Py Se 
2017 年 05 月 底 曾 对 相机 的 安装 位 置 和 伺服 系统 进行 调整 ， 像 场 发 生 了 变化 . 因此 , 使 


TER 


护 前 后 的 相机 姿态 . 


数 . 图 3 显示 了 其 用 户 界面 , 在 左上 角 列 出 了 茶 选 定 恒星 
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可 还 原 大 气 层 的 光影 效果 ， 并 冯 加 赤道 式 和 地 平 式 举 标 网 格 , 满足 本 研究 的 需 
经 度 :100°6'32.07”)、 海 拔 (4800 m)、 
并 查找 某 恒星 的 地 平 坐标 值 等 参 


以 及 指定 的 时 间 ， 即 可 


O 195164 - HD 29453 


reae = 
:©@ + ss 


Ka 


件 下 的 地 平 坐 


妈 3 Stellarium 虚拟 天 文 馆 的 / 


为 了 人 


WN. 


条 得 的 4 组 全 


i 于 显示 和 打印 ,已 经 对 星 图 做 了 色调 反 转 处 


天 彩色 图 像 . 其 


E11 帧 ， 


HFM 


全 


， 夜 间 两 组 共 两 帧 , 均 


量 像 场 . 数据 的 基本 参数 列 于 表 1. 于 


28:08 UTC*08:00 


ome AADA "Gb? 


TE 


Fig. 3 User interface of the Stellarium Virtual Observatory. 


3 测量 像 场 的 几何 参数 


该 全 天 相机 指向 天 顶 拍摄 ， 获得 近似 等 角 投影 的 全 天 空 半球 像 ， 天 空 物 点 的 旨 
距 和 极 角 . 在 像 场 的 边缘 ,部 分 天 际 线 被 山 益 、 监 测 设备 等 物 谈 挡 ， 我 们 采用 
中 心 、 半 径 等 几何 参数 


3.1 选取 监测 数据 


无 名 山 全 天 相机 于 2016-10-05 日 和 2017-06-05 HIHET T Et 
集 了 285 和 352 帧 图 片 . 经 观察 ， 在 
山村 和 镜头 的 边框 均 被 照 亮 ,， 太 阳光 


于 对 视 场 边缘 进行 定位 ,例如 医 


3.2 测量 方法 和 结果 


中 的 均值 . 在 后 


对 
Æ, MWK, HF 


如 前 所 述 , 在 Dl-d 和 D3-d 的 


位 对 半径 、 


业 对 


ba 


将 像 场 的 边缘 点 进行 椭圆 拟 合 ， 
续 计算 


1. 我 们 选取 了 两 组 


将 基于 这 些 参数 对 图 


青天 太阳 或 满 


HARKE 


5 Ft 


共计 11 WEA, 4 


D1-d 和 D3-d 组 数据 的 测量 结果 列 于 表 2. 长 
长 度 的 均值 即 为 


圆 半径 . 


观测 日 


期 之 间 ， 


Huby XY 坐标 分 别 进行 归 一 化 处 理 , 可 得 
0.7590, 0.5048] 和 [0.4295, 0.7611, 0.5044], 


卖 的 监测 ， 使 用 两 种 时 间 和 空 
起 /落下 的 时 刻 ， 入 射 光线 


色 椭 圆 是 拟 合 的 


用 这 两 个 月 的 数据 比较 维 


度 和 方位 角 对 应 着 像 点 的 极 
标记 边缘 点 的 方法 测量 圆 形 像 场 的 


x 间 分 辨 率 分 别 采 


的 天 顶 角 约 为 90"， 视 场 边 缘 的 
也 没有 直射 镜头 而 产生 炫 光 ,， RRR ARR. 而 像 场 外 基本 黑暗 ， 便 
分 别 编号 为 D1-d 和 D3-d. 


长 、 中 心 坐标 、 偏 心率 等 几何 参数 ， 并 计算 它们 在 每 组 图 像 
像 进行 分 析 . 在 图 1 中 , A 


象 场 边缘 . 


半 轴 和 短 半 轴 


的 差异 


曾 对 全 天 相机 进行 了 深入 维 


MARA 


的 千 分 之 


护 . 若 以 图 
D1-d 和 D3-d 的 上 述 三 量 的 归 


Z, 像 场 可 近似 为 区 


bau 


片 Y 方 向 的 分 辩 率 为 单 


一 化 数值 分 别 为 [0.4289， 


二 者 的 差异 分 别 为 [0.14%, 0.28%, -0.07%]. F 
像 场 基本 几何 参数 的 改变 是 微小 的 , 至 少 在 云 量 分 析 中 可 以 忽略 . 


日 此 可 知 ,常规 的 维护 作 


1 http://stellarium.org/ 
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4 测量 天 顶 位 置 


4.1 测量 原理 
鱼 眼 镜头 成 像 最 显著 的 光学 畸变 就 是 圆 形 像 场 中 的 桶 形 像 差 . 假设 鱼 眼镜 头 的 桶 形 像 差 是 轴 对 称 的 ， 并 姑 
是 忽略 像 场 中 心 与 天 项 的 偏差 那么 , 投影 中 桶 形 像 差 只 对 像 点 的 极 距 有 非 线 性 的 影响 ， 并 不 改变 极 角 , 任意 两 
个 恒星 像 点 对 天 顶 投 影 点 形成 的 张 角 是 不 受 桶 形 像 差 影 响 的 ， 应当 与 天 球 上 的 实际 张 角 相 同 . 依据 这 个 原理 ， 
首先 选取 若干 枚 恒星 , 在 图 像 上 测量 其 坐标 ， 同 时 ， 使 用 虚拟 天 文 馆 得 到 其 地 平 经 纬度 . 然后 在 天 顶 区 域 构建 
一 个 网 格 , 逐 格 点 计算 各 恒星 像 点 对 格 点 所 形成 的 张 角 ,并 与 虚拟 值 进行 比较 . 当 实 测 张 角 与 虚拟 张 角 最 为 接 
近 ， 即 二 阶 曲率 为 正 时 ， 则 可 判定 该 格 点 对 应 着 天 项 
在 选取 恒星 时 ， 当 本 着 以 下 三 个 原则 : (1) 低 星 等 . 该 全 天 相机 使 用 的 是 商业 CCD 而 不 是 星 敏 感 器 ， 只 对 
低 星 等 恒星 才 和 全 EME MEM. OTE. 相 比 于 天 文 望 远 镜 ， 全 天 相机 的 像 元 分 辨 率 是 相当 低 的 ， 定 位 精度 约 为 4， 
不 能 有 效 区 分 星团 、 双星 、 或 数 密度 大 的 天 区 . 挑选 狐 立 、 与 背景 有 较 高 反差 的 恒星 ,便于 识别 . (3) 中 等 纬度 . 高 
地 平 纬度 的 恒星 与 天 项 的 角 距 离 较 近 , 在 计算 极 角 时 对 测量 误差 较为 敏感 ; 而 低地 平 纬 度 恒星 光 程 较 长 且 成 像 
受 镜头 像 差 影响 较 大 ， 二 者 均 不 宜 选 作 本 实验 的 样本 . 
基于 上 述 分 析 ， Sie eal 0 大 (~200) 三 个 容量 的 恒星 样本 ， 分 别 测量 了 天 项 位 置 . 下 面 ， 
先 以 小 容量 样本 介绍 计算 步骤 ,然后 介绍 三 个 恒星 样本 的 计算 结果 . 


4.2 选取 恒星 样本 

2016-10-20 日 无 名 山 气象 条 件 较为 理想 ,全 天 相机 工作 正常 ， 获 得 了 时 序 稳定 、 图 像 清 晰 的 时 间 序 列 . 在 日 
落 (18:39”) 与 月 亮 升 起 (22:41) 之 间 ”， 天 空中 没有 强 光源 干扰 ， 名 天 区 的 星空 = 均 较 为 清晰 , 满足 本 研究 的 需要 . 经 
挑选 , 采用 于 21:21:00 获取 的 图 片 (图 4 左 ) 作 为 本 研究 的 样本 ， 记 为 D2-n( 表 1). 

按照 前 文 所 述 的 原则 , 选取 了 9 枚 恒星 , 它们 的 地 平 纬度 在 30 一 60° 左 右 ， 地 平 经 度 分 布 较为 分 散 ， 记 为 样 
本 S2-9. 采用 二 维 高 斯 拟 合 测量 其 坐标 ， 以 下 称 为 实测 星 场 . 图 5 中 ,在 局 部 放大 的 反 转 图 像 上 , 方 片 标记 了 样 
本 像 点 的 位 置 . 同时 , 在 Stellarium 虚拟 天 象 软件 中 ， 获 得 恒星 样本 在 D2-n 拍照 参数 下 的 地 平 坐标 值 ， 以 下 称 之 
为 虚拟 星 场 . 表 3 列 出 了 样本 S2-9 在 图 像 正 交 坐 标 系 中 的 实测 坐标 、HIP (Hipparcos catalogue 中 ) 编 号、 及 其 地 
PHAR. 图 6 在 地 平 坐标 中 绘制 了 样本 的 位 置 . 中 心 十 字 为 天 顶 ; 箭头 标记 正 北方 向 ; 刻度 标记 地 平 纬度， 单 
位 为 度 . 


”北京 时 间 ， 下 同 . 
:是 时 间 序 列 中 日 月 消光 的 时 刻 , 不 是 日 月 到 达 绝 对 地 平 的 时 刻 . 
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图 4 (A) 样本 D2-n, 2016-10-20 H 21:21:00 全 天 相机 获取 的 夜间 图 像 . ( 右 ) 旋 转 后 的 实测 星 图 . 


X /Pixel 


Fig. 4 (left:) Sample d2-n, the night sky captured by the all-sky camera on 21:21:00 of 2016-10-20. (right:) Rotated star map. 
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图 5 恒星 样本 S2-9 像 点 的 分 布 ， 


X /pixel 


方 片 标记 


Fig.5 Distribution of stellar sample S2-9’s image points (marked with spades). 
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6 在 地 平 坐标 系 中 绘制 的 恒星 样本 S2-9( 星 号 ) 的 分 布 


Fig.6 Star sample S2-9 (marked with asterisks) as plotted in a horizon coordinate system. 


4.3 计算 天 顶 位 置 


从 时 间 序 列 中 初步 确定 天 顶 位 于 HIP109056 和 HIP107975 Fy! 
颗 恒星 的 坐标 ， 记 于 表 3 的 末 两 行 . 
Z, 在 图 5 FH H Fr Sk J AD AR a FEB pid. 
在 区 域 Z 构建 一 个 网 格 ， 以 对 天 顶 进行 精确 的 定位 . 在 其 所 有 格 点 中 , 样本 恒星 对 格 点 张 角 的 实测 值 与 虚 


点 附近 . 在 虚拟 和 实测 星 场 中 分 别 测量 这 两 


以 它们 的 中 点 为 中 心 选 取 一 个 边 长 约 为 4? 的 方形 区 域 ， 记 为 天 项 候选 区 域 


拟 值 最 为 接近 的 , 就 是 天 顶点 . 


首先 ， 对 区 域 乙 中 的 每 一 个 网 格 点 (x, y), 计算 第 i1 和 i-1 颗 恒星 像 点 对 


KAG y). 其 次 , 根据 各 恒星 方位 角 的 虚拟 值 ， 1 颗 恒星 对 真实 天 项 所 形成 的 张 角 64. 
N 枚 样本 恒星 张 角 数列 的 虚拟 值 和 实测 值 的 综合 差 
6(x,y) = stddev({0 (x m — 0;}), 1<i<N. 


图 7 绘制 了 某 次 计算 中 网 格 内 5(zx,y) 的 坑 状 分 布 ,X 和 立轴 的 坐标 单位 为 0.02 像素 , 纵 轴 的 单 
二 阶 曲率 为 正 ， 这 说 明天 项 确实 在 乙 中 . 该 坑 的 最 低 点 具有 最 小 的 6(x,y) 值 ， 它 就 是 图 像 的 天 顶 位 置 


其 所 形成 的 


然后 ， 比 较 


WARE. 5(x,y) 的 


考虑 到 ,在 对 6(x,y) 的 数值 计算 中 ,各 恒星 的 排序 会 影响 计算 的 结果 . 这 一 方面 是 因为 张 角 对 中 心 点 位 置 的 


敏感 度 与 张 角 的 大 小 成 正比 ， 另 一 方面 是 因为 对 恒星 像 点 的 测量 有 和 角 分 级 (像素 尺度 ) 的 误差 . 为 了 减 小 某 特 定 
者 果 的 重心 I 在 图 8 所 示 的 计算 结果 中 , REITT 1024 次 计算 ,每 次 使 


排序 中 的 随机 误差 , 将 多 次 计算 牧 


用 随机 函数 对 样本 S2-9 中 的 恒星 进行 随机 排序 ， 共 得 到 1024 个 天 顶 候选 位 置 ， 它们 分 布 在 一 个 4x4 像素 的 矩形 
区 域内 ,候选 天 顶 的 平均 56(x,y) 为 0.39°. UEP EE IB 的 重心 , 格 点 间距 为 0.02 像素 ， 坐 标 单位 为 


BA. ERR 1.5 倍 标准 差 以 外 的 离 群 值 并 取 平 均值 后 ， 得 到 天 项 位 置 : (1294.87, 853.50) 像 素 . 使 用 不 同 的 网 格 
密度 、 离 群 值 剔 除 标准 nc、 计算 次 数 ， 对 样本 S2-9 进行 了 多 次 计算 ， 所 得 天 项 坐标 上 只 有 


像素 级 的 差异 ， 如 表 4 


所 列 . 在 图 5 E, 小 十 字 标 记 从 天 项 候选 区 域 ZX 灰色 方块 ) 中 遵 选 出 来 的 天 项 ; 三 角 号 标记 全 天 相机 圆 形 像 场 的 
中 心 ; 背景 是 图 像 D2-n 的 相应 区 域 ; 坐标 单位 为 像素 . 


4.4 使 用 较 大 容量 样本 的 计算 结果 
在 D2-n 星 图 上 男 选 取 了 一 个 含有 204 枚 恒星 的 样本 S2-204， 并 重新 计算 了 天 顶 坐 标 , 结果 列 于 表 4. 


同 参数 的 多 次 计算 中 , 散 点 均 形成 了 聚合 的 团 禾 ee 只 有 亚 像素 级 的 差异 , 但 是 


在 使 用 不 


天 顶 在 X 


立方 向 与 $2-9 的 计算 结果 均 有 约 3 像素 的 位 移 . We 样本 造成 的 ,对 恒星 像 点 位 置 测 量 的 随机 误差 在 高 
纬度 天 区 得 到 了 放大 并 进而 影响 到 天 顶 的 定位 ,与 计算 参数 的 设置 基 KER. 因为 这 两 个 恒星 样本 的 容量 有 数 


量 级 的 差异 , 我们 认为 S2-204 WA 


直 果 更 可 靠 . 在 表 4 第 10 行 所 列 的 计算 中 ,选用 了 较 小 的 格 点 间距 和 no, EK 
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有 为 天 顶 坐 标 : (1298.03+0.34, 856.60+0.38). 以 图 像 Y 了 方向 的 长 度 为 单位 ， 归 一 化 天 顶 坐 标 为 (0.75118+0.00020， 
0.49572+ 0.00022). 


图 7 AWK Z 网 格 内 6(%,y) 的 坑 状 分 布 
Fig.7 Pit-shaped distribution of 5(x,y) in the Zenith region Z grid. 


++ 


++ 


++ 
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X /pixel 


图 8 使 用 样本 S2-9 计算 所 得 的 候选 天 顶点 


Fig. 8 Distribution of computed candidate zenith, with its centroid marked with a spade. 


4.4 相机 的 竖 直 姿态 偏差 

测量 相机 的 光 轴 与 天 项 的 偏差 .是 本 研究 所 要 解决 的 重要 实际 问题 . 首先 , BEA Z 区 域 的 HIP109056 和 
HIP107975 作为 参考 点 , 依据 二 者 星 点 的 角 距 离 和 像 点 的 距离 , 求 得 相机 在 天 项 区 域 的 平均 空间 分 辨 率 为 0.1145° 
像素 .由 此 可 得 ,图像 的 天 顶 位 置 (1298.03, 856.60) 与 像 场 中 心 (1311.48, 872.22) 的 角 距 离 为 2.36°. 

针对 2017 年 05 月 底 相机 姿态 被 调整 的 情况 , 我 们 使 用 D3-n 星 图 男 采集 了 一 个 含有 40 枚 恒星 的 样本 S3-40， 
代表 性 结果 列 于 表 4 第 11 行 . 天 顶 位置 与 像 场 中 心 的 角 距 离 达 7.64°, 该 次 维护 作业 明显 加 大 了 坚 直 姿态 偏差 . 


5 测量 旋转 角 


全 天 相机 的 旋转 角 w， 即 图 像 纵 轴 与 测 者 子午 线 正 北方 向 的 夹 角 ,可 以 根据 恒星 样本 方位 角 的 虚拟 值 和 实 
测 值 的 差 值 求 得 . 


y = mean({ gi — qi}), 0<i<N. 


202006.00221v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


李 小 波 等 : 基于 虚拟 天 文 馆 对 全 天 相机 姿态 的 精确 测 : 


Hi 


Est, Qi 和 中 DALE SE A 拟 星 场 中 恒星 i 的 方位 角 . S2-9 和 S2-204 的 旋转 角 分 别 为 152.74+2.55° 
和 152.7141.58°. 图 9 绘制 了 S2-204 的 { pi- ph. AS 和 6 中 的 实 线 箭 头 分 别 标记 了 实测 与 虚拟 星 图 的 正 北方 
向 ,舌头 与 Y 轴 的 夹 角 是 图 像 的 旋转 角 . 将 原始 图 像 旋转 y 后 ， 二 者 重合 . 图 4( 左 ) 和 ( 右 ) 分 别 是 旋转 前 后 的 样本 
D2-n. 小 十 字 标 记 天 顶 ， 坐 标 单位 为 像素 . 在 图 4( 右 ) 中 ,图 像 纵 轴 为 正 北方 向 . 这 样 就 制 成 了 符合 行业 规范 的 图 
片 和 动画 . 

在 对 S3-40 样本 的 计算 结果 中 ( 表 4)， 其 旋转 角 为 185.78°， 与 D2-n 的 测量 结果 有 较 大 差异 . 在 近年 的 使 用 中 ， 


每 次 开 箱 维护 都 会 改变 旋转 角 . 因此 , 在 每 次 维护 之 后 都 应 当 重 新 测量 旋转 角 并 反馈 给 数据 处 理 程序 . 
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图 9 用 样本 S2-204 计算 图 像 D2-n 的 旋转 角 


Fig. 9 Calculation of night image D2-n's rotation angle with star samples S2-204. 


结论 与 讨论 


结论 

我 们 使 用 无 名 山 全 天 相机 于 2016 年 10 月 和 2017 年 06 月 获得 的 4 组 全 天 图 像 , 采集 了 容量 分 别 为 9、40、 
204 的 恒星 样本 ， 从 虚拟 天 文 馆 获 得 了 恒星 的 地 平 坐 标 , 通过 像 场 对 姿态 参数 进行 了 精确 的 测量 . 主要 结果 为 : 
O 开发 了 一 种 基于 虚拟 天 文 馆 计算 恒星 坐标 功能 的 测量 方法 ,使 姿态 测量 建立 在 容量 较 大 的 参考 恒星 上 ,为 全 
天 相机 图 像 的 深入 处 理 分 析 黄 定 了 基础 。 特别 是 , 采用 网 格 法 对 天 顶 的 定位 达到 了 亚 像素 的 精度 ， 有 望 满足 陨 星 
和 航空 器 等 高 定位 精度 监测 的 需求 . 该 方法 对 相机 的 视 场 、 像 差 、 可 见 天 空 角 面积 等 参数 没有 依赖 性 ， 具 有 较 
强 的 可 移植 性 。 
(2) 提高 相机 姿态 精度 的 首要 任务 是 减 小 光 轴 的 天 顶 角 ,而 对 相机 的 维护 作业 可 能 会 造成 光 轴 与 天 顶 有 不 可 忽 
略 的 偏差 . 按 2016-10-20 日 全 天 相机 的 姿态 , 像 场 的 天 项 坐标 为: (0.75118+0.00020, 0.49572+0.00022)， 与 像 场 中 
心 的 夹 角 为 2.36°; 而 在 2017-06-05 的 数据 中 ， 二 者 的 夹 角 达到 7.6°. 
(3) 常规 维护 作业 对 像 场 几 何 参 数 的 改变 其 微 (<0.28%)， 在 一 般 图 像 分 析 中 可 以 忽略 . 无 名 山 全 天 相机 的 像 场 呈 
正 圆 形 . 以 图 像 Y 了 方向 的 长 度 为 单位 其 半径 和 中 心 坐标 分 别 为 0.4289 和 (0.7590、0.5048). 
(4) 测量 了 图 像 的 旋转 角 并 在 可 视 化 产品 中 按 制图 惯例 进行 了 旋转 . 每 次 维护 相机 都 可 能 改变 旋转 角 ， 并 需要 重 
新 测量 


6.2 讨论 


本 测量 的 主要 误差 来 源 是 : (1) 单 帧 图 像 中 恒星 像 点 的 定位 误差 。 夜 间 长 时 间 曝 光 (~30s) 和 视 宁 度 的 时 空 不 
均匀 性 等 均 可 造成 偶然 误差 . (2) 光学 系统 的 桶 形 像 差 。 由 于 光 轴 倾斜 ， 光 学 投影 的 径 向 非 线性 被 传导 到 切 向 ， 
在 测量 结果 中 形成 了 系统 误差 . 在 进一步 的 工作 中 ， 可 以 采用 升级 硬件 以 缩短 曝光 时 间 、 综 合 多 帧 测量 结 
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改进 测 星 方法 等 手段 抑制 偶然 误差 ， 通 过 测量 镜头 径 向 畸变 参数 ， 在 预 处 理 中 将 图 像 矫正 为 等 角 投 影 ， 以 抑制 
系统 误差 . 
影响 全 天 相机 像 场 的 首要 姿态 偏差 是 其 光 轴 与 天 顶 的 夹 角 . 如 果 光 轴 不 坚 直 ,不 仅 对 视 场 和 投影 的 轴 对 称 
性 都 会 造成 系统 性 的 影响 ,而 且 每 次 在 维护 后 都 需要 重新 测量 姿态 ,如 果 能 在 设计 和 施工 中 大 幅度 提高 坚 直 方 
向 的 姿态 精度 和 稳定 性 ,在 野外 使 用 中 光 轴 的 天 顶 角 一 直 符 合 监测 任务 的 精度 要 求 ， 那 么 , 只 需要 测量 新 的 旋转 
角 就 可 以 了 . 
在 天 文 观测 中 ， 一 个 潜在 的 误差 来 源 是 观测 设备 的 时 钟 误差 .由 于 天 球 旋转 的 速度 是 360°d-!1 即 15”s-1， 
时 钟 误差 0.1 s 就 会 造成 角 秒 级 的 定位 误差 . 但 是 ， 对 于 无 名 山 全 天 相机 而 言 ， 按 照 天 顶 区 域 0.0572° 的 像 元 分 
WER (5184x3456) ， 时 钟 误差 达到 13.7 s 方 能 造成 像素 级 的 定位 误差 .无 名 山 全 天 相机 使 用 互联 网 时 钟 对 工控 
机 内 置 时 钟 进行 校准 ， 其 精度 完全 可 以 满足 恒星 定位 的 需要 . 
天 文 导航 技术 的 发 展 , 为 本 研究 的 进一步 深入 提供 了 很 多 可 借鉴 的 方法 . 本 研究 与 天 文 导 航 的 研究 具有 一 
定 的 异同 . 一 方面 , 天文 导 航 的 测量 目的 是 定位 ， 要 求 准 实时 性 ; 而 本 研究 的 测量 目标 是 地 基 设 备 的 姿态 ,对 时 
效 性 要 求 不 高 ， 但 期 望 在 类 似 设备 中 具有 较 强 的 可 移植 性 . 另 一 方面 ,二 者 都 主要 解决 实测 星 场 的 获得 、 矫 正和 
比较 计算 . 全 天 相机 充足 的 恒星 样本 容量 , 提高 了 结果 的 可 靠 性 ， 并 为 深入 分 析 提 供 了 可 能 . 如 果 在 后 续 研 究 中 
需要 采集 大 容量 的 恒星 样本 ,可 以 从 天 文 导航 的 研究 中 借鉴 成 熟 的 从 星 图 中 批量 获取 和 比较 恒星 坐标 的 方法 
273] 并 开发 “从 虚拟 天 文 馆 中 自动 输出 坐标 矩阵 的 接口 软件 . RIERS AT ee 
法 ,开展 对 全 天 相机 光学 系统 的 分 析 和 矫正 ， 并 提高 对 数据 格式 、 图 像 瑕 辣 和 气象 条 件 的 适应 能 力 ， 以 期 把 本 而 
究 移植 到 类 似 设备 . 


致谢 无 名 山 观测 站 在 次 仁 汪 堆 、 洛 绒 次 称 、 宋 其 武 等 同志 的 辛苦 劳作 下 得 以 稳定 运行 ,在 此 深 表 谢 忱 


4 Stellarium 软件 执行 的 是 GNU General Public License, 用 户 可 以 修改 程序 . 
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表 1 本 研究 使 用 的 全 天 相机 监测 数据 
Table 1 Data of all-sky camera used in this study 


期 时 间 。 分 辩 率 间隔 时 长 帧 数 天 气 测量 对 象 编号 
2016 年 10 月 05 06:29:11-19:52:38 2592x1728 5.05 min 7 多 云 间 晴 视 场 边界 Dl-d 
2016 年 10 月 20 日 21:21:00 2592x1728 5.05 min 1 晴 相机 姿态 D2-n 
2017 年 06 月 05 02:20:55-22:39:22 5184x3456 4.09 min 4 晴 转 多 云 视 场 边界 D3-d 
2017 年 06 月 05 日 03:54:58 5184x3456 4.09 min 1 晴 相机 姿态 D3-n 

(a) 对 含有 多 帧 图 像 的 数据 组 ， 此 处 列 出 起 止 帧 的 时 刻 , 下 同 


O) 全 天 相机 在 该 日 采集 时 间 序 列 的 时 间 间 隔 . 


表 2 全 天 相机 像 场 的 实测 几何 参数 


Table2 Geometric parameters of all-sky camera's image field 


数据 编号 分 辨 率 长 半 轴 “ 短 半 轴 圆 半 径 Pate Aas 中 心 纵 坐标 
Dl-d 2592x1728 741.25 741.02 741.14 1311.48 872.22 
D3-d 5184x3456 1485.63 1483.18 1484.40 2630.20 1743.16 
(a) 长 度 单位 均 为 像素 . 
表 3 恒星 样本 S2-9 的 实测 和 虚拟 参数 
Table 3 Measured and virtual parameters of star sample S2-9 
序列 号 坐标” 立 坐 标 HIP 编号 地 平 经 度 地 平 纬度 
0 1641 1236 113368 165°07'45.6" 29°40'22.6" 
1 1797 1012 3419 134°26'36.5” 28°21'44.1"" 
2 1264 1285 100345 211°27'04.7" 40°16'36.3" 
3 1293 519 746 27°02'28.4" 52°32'56.7" 
4 1679 569 9884 80°07'24.3” 34°34'05.2"" 
5 1588 797 1067 106°37'54.3” 56°00'25.4" 
6 1110 1114 97649 242°25'22.5" 53°55'30.6" 
7 1111 766 102098 322°08'55.6" 67°45'51.9"" 
8 919 855 91262 296°59'40.3" 47°37'09.5" 
—° 1281 856 107975 260°11'39.3” 88°41'05.5" 
= 1303 844 109056 91°43'00.0"" 88°25'55.4" 


(a) xX. Y 是 图 像 D2-n 上 恒星 像 点 的 正 交 坐标 , 单位 为 像素 . 
(b) 样本 恒星 在 图 像 D2-n 拍摄 时 刻 的 地 平 经 纬度 ， 由 虚拟 天 文 馆 计算 所 得 . 
(c) HIP107975 和 HIP109056 仅 用 于 圈定 天 项 的 位 置 范 围 , 在 计算 天 项 时 并 不 作为 参考 恒星 . 
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表 4 天 顶 和 旋转 角 的 计算 结果 


Table 4 Results of the computation of zenith and rotation angle 


恒星 样本 计算 参数 计算 结果 

编号 星 图 容量 。 格 点 间距 。 n? HARMO RWM 天 项 纵 坐 标 旋转 
0.50 2.0 512 1294.91 853.43 152.75 
0.20 2.0 512 1294.84 853.50 152.74 
S2-9 D2-n 9 0.10 20 512 1294.82 853.41 152.74 
0.02 1s 1024 1294.87 853.50 152.74 
0.01 20 32 1294.99 853.61 152.74 
0.50 1.5 1024 1298.04 856.64 152.71 
0.10 15 256 1297.93 856.47 152.71 
S2-204 D2-n 204 0.05 15 512 1298.02 856.53 152.71 
0.02 5.0 1024 1297.65 855.66 152.71 
0.02 1.5 1024 1298.03 856.60 152.71 
$3-40° D3-n 40 0.05 3.5 256 2561.85 1627.99 185.78 


(a) 格 点 间距 和 天 项 坐标 的 单位 均 为 像素 ,旋转 角 的 单位 为 度 . 

(b) n 是 剔除 离 群 值 的 标准 . 距离 散 点 重心 大 于 na 的 点 被 视 为 离 群 点 ，a 是 标准 差 . 
O 每 次 计算 中 , 对 样本 恒星 进行 随机 排序 ,结果 是 剔除 离 群 点 之 后 的 均值 . 

(d) 仅 列 出 一 代表 性 结果 . 
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